Effects of spatial non-uniformity on laser dynamics. 
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Semiclassical equations of lasing dynamics are re-derived for a lasing medium in a cavity with 
a spatially non-uniform dielectric constant. The non-uniformity causes a radiative coupling be- 
tween modes of the empty cavity, which results in a renormalization of self- and cross-saturation 
coefficients. Possible manifestations of these effects in random lasers are discussed. 
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Introduction Random lasers, in which optical feed- 
back is provided by scattering of light due to spatial inho- 
mogeneity of the medium rather than by well defined mir- 
rors, has recently attracted a great deal of attention (HQ 
In the case of weak scattering, when the propagation of 
light can be described within diffusion approximation, 
the nature of lasing in such systems has been well under- 
stood starting with a pioneering work by LetokhovQ fol- 
lowed by a large volume of subsequent experimental and 
theoretical studies. The case of strong scattering, how- 
ever, when light can be at the verge of Anderson localiza- 
tion, remains much more controversial. Experimental re- 
sults of Ref . |4{ and consecutive works with strongly scat- 
tering systems (see recent reviews in RefQ and Ref.Q) 
led to an assumption that lasing observed in those ex- 
periment is due to formation of pre-localized, if not com- 
pletely localized states of light, which play a role of lasing 
cavities 5], and provide coherent resonant optical feed- 
back as oppose to non-resonant feedback affecting only 
intensity of light in the diffusion case. The presence of 
narrow multiple lasing peaksQ as well as Poisson statis- 
tics of emitted radiation^ were considered as evidences 
in the favor of this interpretation of these experiments. 
However, it was shown in Ref.Q that the non-resonant 
feedback can also result in lasing with multiple narrow 
peaks. Moreover, the authors of Ref.Q demonstrated 
that the Poisson statistics also cannot be considered as 
an exclusive attribute of lasing with the resonant feed- 
back. 

In this situation, the recent r esults of Ref© assume 
a particular significance. In these experiments a multi- 
peak lasing was observed in PMMA sheets containing a 
rhodamine dye as an active material and titanium diox- 
ide microparticles as scatterers. This system is char- 
acterized by a strong inhomogeneous broadening of the 
lasing transition, and most of the lasing peaks are sep- 
arated in the frequency domain by a homogeneous line 
width of the lasing transition, j a . This is naturally ex- 
plained by the competition of modes " feeding" from the 
same population inversion and spectral hole burning in 
inhomogeneously broadened systems^- However, above 
certain value of the pumping intensity, there were ob- 
served two lasing peaks coexisting within the homoge- 
neous line width, j a and having synchronized temporal 



behavior . This observation is indicative of the genuine 
two-mode lasing, which can occur in regular cavity lasers, 
when the mode competition is weakened by spatial hole 
burning 10] . Such a behavior, however, cannot take place 
in the the case of the non-resonant feedback, because, in 
diffusive systems lasing only occur at the frequency of 
an atomic transition 3] (multiple peaks in Ref. [7| are due 
to inhomogeneous broadening of the transition used to 
generate emission and do not signify a truly multi-mode 
behavior). 

Thus, as of today, the results of Ref fcj provide the most 
convincing evidence of the resonant feedback in random 
lasers. It is important, therefore, to achieve a clear un- 
derstanding of the specifics of non-linear mode interac- 
tion in such systems. However, since the experiments of 
Ref|3 deal with just a single realization, the randomness, 
by itself, is not important here. What is important is 
the spatial inhomogeneity of the quasi-cavity, supporting 
the modes of interest. The main objective of this paper, 
therefore, is to study how this spatial inhomogeneity af- 
fects lasing threshold and spatial hole burning. Our con- 
sideration, however, is not constrained by random lasers, 
and can be applied to any type of lasers with spatially 
non-uniform cavities. Currently, there is a tremendous 
interest to lasing in the systems with a modulated di- 
electric constant, for instance, photonic crystals. The 
result presented here are relevant for these systems as 
well. Moreover, current technologies allow for engineer- 
ing structures with virtually arbitrary spatial profile of 
the dielectric function. The results of this paper can be 
used to manipulate properties of lasers by using spatial 
dependence of the cavity dielectric function as a new de- 
sign parameter. 

A general multi-mode theory of lasing in systems with 
an arbitrarily inhomogeneous dielectric constant, e(r), 
presented here is an extension of semi-classical Lamb 
theory^ for the media whose dielectric constant is in- 
homogeneous in the direction of propagation of the laser 
beam (inhomogeneity in the perpendicular directions re- 
sults in waveguiding effects, which are well studied in 
laser physics (see, for instance, 01) • This inhomogeneity 
modifies the orthonormalization condition for the eigen 
modes of the cavity, making the standard inner product 
of the modes belonging to different eigen frequencies dif- 
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ferent from zero. The main effect resulting from this non- 
orthogonality is a new type of linear coupling between 
normal modes of the empty cavity, which is mediated by 
the polarization of the active medium. 

The non-orthogonality of eigen modes due to the in- 
homogeneity of e(r) should not be confused with non- 
orthogonality of Fox-Li modes of uniform but leaky 
(open) cavities, which arises due to non-hermitian nature 
of the respective eigen value problem and does not result 
in any additional linear coupling between the modes [To| . 
The only consequence of the non-hermitian nature of 
such cavities is the presence of an additional factor in 
the linear susceptibility of the active medium, which was 
carefully studied in the past and shown to be responsible 
for the excess noise in unstable cavities[l2j. 

Multi-mode laser equations for an inhomogeneous 
medium We consider an ideal cavity specified by an in- 
homogeneous dielectric function e(r) and some boundary 
conditions. The cavity is filled with an active medium 
characterized by its polarization P(r). Let us assume 
that we know the full system of eigen modes, /fe(r), and 
respective eigen frequencies ui k of such a cavity in the 
absence of the polarization. These modes can be used to 
present electric field, E, and polarization P in the form 
of their linear combinations: E = ^ fc Ek(t)fk(r), P = 
J\ Pk (t)fk (r), where we assume that only s-polarizcd 
modes couple to the active medium, and ignore the vector 
nature of the field and the polarization. The orthonor- 
malization condition for these modes involves inhomoge- 
neous dielectric function e(r)0:J e(r)/^ (r)/fc 2 (r)cZr = 
Ski k 2 ! which means that the wave functions f k (r) them- 
selves are neither normalized nor orthogonal. As a result, 
the dynamic equations for the amplitudes, Ek, takes the 
following form 



E k (t) + (aJk-i-fk) 2 Ek(t) 



-AnJ2v kkl P kl (t), (1) 

fei 



where we introduced cavity losses, characterized by phe- 
nomenological parameters jk- The main peculiarity of 
Eq. Q is the presence of the linear coupling between dif- 
ferent polarization amplitudes, P k , characterized by non- 
diagonal elements of the matrix 



(2) 



V kkl = / f* kl (r)f k2 (r)dT 



The presence of such a coupling is the main difference 
between homogeneous and inhomogeneous cavities. The 
magnitude of coupling parameters V kkl depends on the 
spatial profile of the dielectric constant, and can be tai- 
lored to enhance (or diminish) the coupling effects. 

Similar equation can be in principle derived for inho- 
mogeneous open cavities as well, where eigen modes f k 
should be replaced by appropriate Fox-Li modes. The 
hermitian orthogonality in this case is replaced by the bi- 
orthogonality, which involves adjoint set of modes. For 



a general case of non-uniform o pen cavity this condition 
was derived, for instance, in RefjlJ. We shall leave, how- 
ever, this topic for future work. 

A gain medium is described within the model of two- 
level atoms, characterized by dephasing rate and 
population relaxation time r, and we use a standard 
density matrix approach in order to derive equations 
for polarization amplitudes P kl and population differ- 
ence AN. The next standard step in the derivation of 
rate equations would be rotating wave and slow ampli- 
tude approximations, which amount to presenting mode 
amplitudes E k and P k as E k (t) = E k (t) exp(—iQ k t), 
P k (t) = H k (t) exp(— ffifct), where E k and life are slowly 
changing amplitudes, and £l k is a frequency of the re- 
spective lasing mode. However, forcing this procedure 
onto equation yields linear oscillating terms of the 
form VkktX^H^) exp[-i(Q k ~ ^fci)] s fei ,which ren- 
der derivation of meaningful rate equations impossible. 
Here 
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(3) 



is a linear susceptibility of the gain medium with \fj,\ being 
dipole matrix element of the lasing transition. Parameter 
ANq represents non-saturated population inversion, and 
characterizes the strength of the pumping. 

The physical origin of this problem is quite clear — in 
the presence of the linear coupling the modes of a pas- 
sive medium are not genuine normal modes of the entire 
system. As a result, an attempt to excite such a mode 
leads to exchange of energy between coupled modes and 
to non-stationary oscillations of the respective intensities. 
The rate equations, therefore, should be derived for the 
normal modes of the entire system, which would include 
cavity and the gain medium. To this end, it is conve- 
nient to transform Eq.JQl in the frequency domain using 
conventionally defined Fourier transformation: 

Sfe [i^k - ilk - S kkl - 2irLu V kkl x (1) (w)] E kl (w|4) 

= 2 7 ro;oE fcl ^fc 1 ^ ) M 

where a tilde on the top of a symbol signifies the Fourier 
transform of the respective quantity, and the polariza- 
tion is separated into a linear and third-order non-linear 
contribution, P^K The former is taken into account in 
Eq.Q by introducing a linear susceptibility an d 
the expression for the latter was derived in a standard 
way from the full system of density matrix equation us- 
ing standard perturbation approach. 
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Anticipating the future use of the rotation wave approx- 
imation applied to genuine normal modes of the system 
(see Eq.© below) I substituted 2lu(uj — uj k + i^ k ) instead 



of i 



[u>k — *7fc) , neglected the non-resonant part of 



the linear susceptibility, and replaced all frequencies uo in 
non-resonant expressions with atomic frequency ujq. The 
latter approximation is justified because we will only con- 
sider the case where frequencies of all participating modes 
lie within a homogeneous line width of the lasing transi- 
tion. Non- linear coupling parameters in Eq.©, Akkik 2 k 3 > 
are defined as 

Akk lk2 k 3 = J e(r)A!(r)/ fcl (r)/ fc2 (r)A;(r)d 3 r (6) 

Lasing threshold and non-linear dynamics of the inten- 
sities In order to illustrate the effects of spatial inhomo- 
geneities on lasing threshold we find eigen frequencies of 
linearized Eq.Q in a two-mode case. Imaginary parts 
of these frequencies, ri j2 are both positive below lasing 
threshold. With increasing pumping, however, one of 
them, Ti, for instance, first changes its sign, and this 
point determines the lasing threshold. If we assume that 
7i *C 72 , a simple expression for the lasing threshold can 
be derived: 



ANS 
AN Qtr = °" 
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(7) 



where N® is a threshold value of ANq in a system with 
a uniform dielectric constant. Two effects of the non- 
uniformity appear in this expression. First, factor Vn, 
which would be equal to unity for a uniform medium, af- 
fects the threshold even in the absence of the linear cou- 
pling between the modes. The value of this parameter 
depends on the spatial profile of the dielectric function; 
with an appropriate choice of the latter one can achieve 
a decrease in the lasing threshold. The second effect re- 
flected in Eq. Q is due to the coupling between the modes 
and results in further decrease of the threshold. 

In order to derive rate equations we have to diago- 
nalize the linear part of Eq.Q. To this end, we will, 
first, neglect the dispersion of the linear susceptibility. 
This approximation is justified if we are only interested 
in dynamics of intensities rather than lasing frequencies, 
and because all the frequencies of interest lie within the 
width of the atomic transition. After that we have to 
solve the eigen- value problem for the remaining matrix, 
which is, however, essentially non-hermitian. Therefore, 
we have to find two adjoint sets of vectors — right (| e,)) 
and left ((ij ), which obey the bi-orthogonality condi- 
tion (efj | ei) = when i =/= j. In order to preserve 
standard expressions for intensities we shall normalize 
our right eigen- vectors using condition (e* | e») = 1 (so 
called power normalizationflfj). As a result the product 
014 = {Si | e,) 7^ 1. In order to eliminate linearly coupled 
terms from Eq.Q we present cavity mode amplitudes, 



Ek, as a linear combination of the right eigen vectors, 
I ei), 

E k {oj) = 27r^T fci [Zi(w)8(u - Qi) + Z,{uj)S(lj + a,)} , 

(8) 

where columns of matrix T k i are formed by the vectors 
| ei). In Eq.JHJ we also introduced a slow changing ampli- 
tude approximation applied to the amplitudes of the true 
normal modes of the system. In the frequency domain, 
this approximation consists of presenting the amplitudes 
as a product of a frequency dependent part and a delta 
function, containing a lasing frequency fij. Matrix T k i 
is not a unitary matrix, and, transformation of Eq.Q 
to the new basis has to rely on matrix Tfo, whose rows 
consist of vectors (£j | . The product of matrices T and 
T is a diagonal matrix with elements a*. The resulting 
equations for slow amplitudes Zi take the form 



.dZi 
' dt 



(fi< - fii) 



nAz t 



(3) -i(n^-n 4 )« 



(9) 

where u tj = J2k,k ± TikV kkl T kl3 j on, and the non-linear 
contribution to polarization, in the new basis, is given by 



N 4 AiVoT 
8^ 3 



i-i{Q,j-Vli+Q, m -Q,i)t 



Hjlm 



j,l,m 



— f2; + f2„ 



UJQ 



(10) 
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(ii) 

where we have neglected frequency dependence of the 
nonlinear coefficients Riji m - Equations |J3J, IjlUI) . and 
(|llfl provide a basis for further analysis of the non-linear 
dynamics of the system under consideration. 

In particular, the rate equations can be obtained in a 
standard way by separating real and imaginary parts of 
Eq.© and neglecting all oscillating terms on its right- 
hand side: 



dh 
dt 



2L 



-r< - - 1 



(12) 



Here Ii is the dimensionless intensity of the i-th mode, I\ 
is its unsaturated amplification rate, fa and 6^ are self- 
and cross-saturation parameters respectively, which are 
expressed in terms of coefficients Riji m of Eci. (|llfl . These 
equations have the same form as standard lasing rate 
equations for a uniform medium, with the only difference 
being that instead of the combination of non-saturated 
gain and loss terms, we have a single parameter Ti, repre- 
senting the imaginary part of the mode's eigen frequency. 
The main effect of the linear coupling is a renormaliza- 
tion of the non-linear coupling coefficients. The most 
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FIG. 1: Dependence of the non-linear coupling parameter C 
on AiVo (curve 1) and deviation of the intensity in a single 
mode regime from a linear dependence on pumping (curve 2) 

important feature of this renormalization, which makes 
it experimentally relevant, is a non-trivial dependence of 
the new coefficients on the intensity of pumping. This 
dependence arises because the coupling is carried by po- 
larization, and, hence, its strength is proportional to the 
unsaturated inversion AA^. In order to illustrate the 
last point we consider an example of only two interacting 
modes. There are two possible regimes of behavior in this 
case: single mode lasing, when the mode competition pre- 
vents the second mode from lasing, and two-mode lasing, 
when the spatial hole burning prevails over the competi- 
tion. The choice between these regimes is determined by 
a coupling parameter C, defined by the ratio of cross- and 
self-saturation parameters C = (610601 ) / (3i 1 800) [Tojj . In 
a uniform medium, this parameter depends solely upon 
spatial distribution of the cavity modes, determined by 
the cavity's geometry. In the situation considered, here, 
this parameter becomes dependent on the pumping in- 
tensity. In order to illustrate the possible character of 
such dependence, we simulated parameter C for a cav- 
ity which consists of two dielectric materials with dif- 
ferent dielectric constants ei, and to.. The curve 1 on 
the figure shows the dependence C(AAo) for two closest 
in frequency modes of such a cavity. The most striking 
feature of this graph is the steep decrease of this coeffi- 
cient with AA*o, which means that even if the modes of 
the empty cavity wouldn't favor the spatial hole burning, 
the increasing with pumping linear interaction between 
the modes modifies their spatial structure in a way which 
is beneficial for the two-mode lasing. Similar behavior of 
C(ANo) was also found for the dielectric constant of the 
shape e(z) = 6q + az 2 or e(z) = eo + Scosz, where z is 
coordinate in the beam propagation direction. This ef- 
fect might explain the two- mode behavior observed in 
random lasersjl"^ . The fact that increased pumping can 



systematically drive C below unity for various configura- 
tions of e(r) makes such effect much more likely to occur 
in a random system than just a coincidental combination 
of various parameters suggested, for instance in Ref . • 
A presence of the linear mode coupling in random lasers 
can be verified directly by observation of dependence of 
I(ANq) in single and multi-mode regimes, which, if ef- 
fects considered here are responsible for the observed 
multi-mode behavior, should deviate from a simple linear 
behavior expected in lasers with the uniform dielectric 
constant (see curve 2 in Fig.l). 

Conclusion We derived non-linear equations describ- 
ing dynamics of lasing modes in a cavity whose dielec- 
tric constant is spatially non-uniform in the direction of 
beam propagation. For a number of spatial profiles of 
e(r) it is shown that the non-uniformity enhances spatial 
hole burning and promote two-mode lasing. This effect 
can explain recent observation of two-mode behavior in 
random lasers. The equations derived in the paper can 
also be used to manipulate properties of lasers through a 
design of spatial profile of the dielectric function. 
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